Self-adjoint wave equations for dynamical perturbations of self-gravitating fields 



O 

o 
o 

(N 
•*-> 
O 

O 



> ■ 
: 
. 

o ■ 
o ■ 

o : 
o ■ 

i ' 



X 



0. Sarbach*, M. Heusler* t , and 0. Brodbeckt 
* Institute for Theoretical Physics, University of Zurich, CH-8057 Zurich, Switzerland 
Time-steps GMbH, CH-8910 Zurich, Switzerland 
(February 7, 2008) 

It is shown that the dynamical evolution of linear perturbations on a static space-time is governed 
by a constrained wave equation for the extrinsic curvature tensor. The spatial part of the wave op- 
erator is manifestly elliptic and self-adjoint. In contrast to metric formulations, the curvature-based 
approach to gravitational perturbation theory generalizes in a natural way to self-gravitating mat- 
ter fields. It is also demonstrated how to obtain symmetric pulsation equations for self-gravitating 
non-Abelian gauge fields, Higgs fields and perfect fluids. For vacuum fluctuations on a vacuum 
space-time, the Regge- Wheeler and Zerilli equations are rederived. 



I. INTRODUCTION 

Perturbation theory finds several interesting applica- 
tions in general relativity. A prominent example is the 
close limit approximation to black hole collisions 0] . Be- 
sides providing estimates for the energy emitted by grav- 
itational radiation, the close limit results also play an im- 
portant role in testing existing numerical codes for full 
non-linear general relativity. Further examples comprise 
the linear stability analysis of neutron stars or black holes 
with matter fields. Recently, a perturbation approach 
has also been used to find new stationary solutions and 
to prove local uniqueness results for hairy black holes and 
self-gravitating solitons §j f| f|. 

For vacuum gravity perturbation theory is highly de- 
veloped, and even for self-gravitating systems important 
general properties of the perturbation equations were es- 
tablished. In particular, it has been shown that a man- 
ifestly hyperbolic, gauge-invariant formulation of gravi- 
tational perturbation theory does exist for static back- 
ground configurations ||, ||. On the other hand, it 
is known for a wide class of matter models that the 
pulsation equations governing spherical perturbations 
of static, spherically symmetric field configurations are 
manifestly symmetric and derivable from the second vari- 
ation of the Komar mass 0. However, only for spe- 
cial systems, such as vacuum gravity and the Einstein- 
Maxwell system for a spherically symmetric background, 
formulations are known for which the perturbations are 
governed by a wave equation which is hyperbolic and 
symmetric. Symmetric formulations are most valuable, 
since they provide the possibility to discuss stability is- 
sues by rigorous means, applying methods from spectral 
theory. 

The commonly adopted procedure to derive the pertur- 
bation equations for spherically symmetric black holes is 
to first use the symmetry of the background in order to 
expand the metric in terms of spherical harmonics. Then, 
the pulsation equations are derived for gauge-invariant 
quantities. For vacuum perturbations of vacuum grav- 
ity, this leads to the Regge- Wheeler equation || in the 
odd-parity sector and to the Zerilli equation || for even- 
parity perturbations. Both equations have the form of 



symmetric wave equations and, as an important conse- 
quence, the stability of the Schwarzschild metric within 
linear perturbation theory can be established. This result 
was extended to the Einstein-Maxwell system by Mon- 
crief |lC| ]. While also starting with an expansion of the 
metric in spherical harmonics, Moncrief uses the ADM 
Hamiltonian to find a symmetric wave equation. In par- 
ticular, he introduces gauge- invariant quantities and at 
the same time, separates the constraint variables from 
the dynamical variables. 

As we have argued in a recent letter [fPH , the metric 
approach fails to yield hyperbolic and symmetric equa- 
tions in the presence of general gravitating matter fields: 
In order to obtain a wave operator appearing already off- 
shell, one has to introduce amplitudes which are adapted 
to the stationarity rather than the spherical symmetry 
of the background. Using ideas introduced by Choquet- 
Bruhat et al. for full, non-linear general relativity ||, 
we show that a hyperbolic, symmetric wave equation 
for the extrinsic curvature can be obtained by natural 
means when the background is static. In contrast to the 
traditional, metric based approach, our curvature based 
approach is manifestly gauge-invariant, hyperbolic and 
symmetric without using the spherical symmetry of the 
background. More importantly, we also show that a nat- 
ural generalization to self-gravitating fields is possible, 
including non-Abelian gauge fields and perfect fluids. 

This work is organized as follows: In Section II, we 
show how to obtain a manifestly hyperbolic and sym- 
metric wave operator for the perturbation of the extrinsic 
curvature on a static background. Then, in Section III, 
we specialize the result on a vacuum background and dis- 
cuss the initial value formulation and the projection onto 
the constraint manifold. The coupling to Yang-Mills- 
Higgs fields and to perfect fluids is discussed in Sections 
IV and V, respectively. In Appendix A, we show how 
to separate the constraint and dynamical variables and 
rederive the Regge- Wheeler and Zerilli equations in quite 
a natural way. Finally, in Appendix B, we recall a useful 
field theoretical formulation of perfect fluids and show 
how this formulation applies to perturbation theory. 
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II. THE CONSTRUCTION OF THE WAVE 
OPERATOR 

The goal is to linearize the equations governing self- 
gravitating matter fields for a static background and to 
bring the perturbation equations into the form of a sym- 
metric wave equation, 



(d t 2 + A) u = 0, 



(1) 



where u describes the perturbed gravitational and matter 
fields, and where A is a (formally) self-adjoint, elliptic op- 
erator containing spatial derivatives up to second order. 
In particular, this implies that equation (|l|) is hyperbolic, 
which guarantees a well-posed initial value formulation. 
Furthermore, as the spatial operator A is self-adjoint, 
methods form spectral theory can be applied and stabil- 
ity issues can be discussed analytically. 

In order to give some motivation for the ideas devel- 
oped below, we construct the wave operator belonging to 
the source-free Maxwell equations in flat space-time. In 
view of a generalization to non-Abelian gauge groups it 
is convenient to formulate Maxwell's equations in terms 
of the electric one-form, E, and the magnetic potential 
one-form, A, where B = *dA. Maxwell's equations con- 
stitute a system of constrained evolution equations, with 
Gauss constraint 

SE = 0, 

and evolution equations 

A = d<P- E, E = SdA. 

Here, the codiffcrcntial operator d> for a p-form u> is de- 
fined by d*u> = (— l) p * d * uj. Note that there is no evo- 
lution equation for the electric scalar potential <p. This 
reflects the gauge freedom of the theory. Usually, a wave 
equation for A is derived by imposing the Lorentz gauge 
condition <j)+d^ A — 0. However, A is not gauge-invariant. 
In order to obtain a hyperbolic wave equation in terms 
of the gauge-invariant quantity E, we differentiate the 
second evolution equation with respect to i, eliminate A 
using the first evolution equation and add the differential 
of the constraint equation. This yields 

E+ (dU + dd^) E = 0, 

which is of the desired form The initial value prob- 
lem may be solved as follows: First, any function <p{t, x) is 
chosen. Next, initial data A t= o and E t= o are given, where 
E t= o is subject to the Gauss constraint. Then, E t= o 
is computed from the second evolution equation, where 
d^E^o = automatically follows. In a next step, E t is 
computed for all times using the symmetric wave equa- 
tion. Finally, the magnetic potential is obtained from 



A t = A t=0 + / (d(/) T - E T ) dr. 



The Gauss constraint and the second evolution equation, 
which was differentiated in order to get the wave opera- 
tor, are both constraint equations for the initial data. It 
is easy to see that they propagate. 

A gauge-invariant wave operator governing linear fluc- 
tuations of self-gravitating fields on a static background 
can be constructed in a similar manner. First, we notice 
that a manifestly symmetric formulation of the pertur- 
bation equations for a vacuum space-time exist: 



i 



(2) 



Here, G^ v is the Einstein tensor, r\ is the volume form on 
M and the tensor G a ^^ v is given in terms of the Riemann 
tensor, R a p^, the Ricci tensor, R@ v = R^a^, and the 
Ricci scalar, R = R^, by 

+ i (2 ff Q V y + 2g^R af3 - 3g ati R 0v - ig^R 01 ") 

+\{g m 9^ + 9^ 9^-9^9^) R- 

Clearly, the operator on the right-hand side (RHS) of 
Eq. (||) is formally self-adjoint with respect to the inner 
product 



(*5 (1 W 2) > = 



g afl g "c 



(3) 



In order to get an evolution equation for the perturbed 
geometry on a static background, our program is the fol- 
lowing: 

1. First, we perform a 3+1 decomposition of Eq. (|^) 
for a static background. The resulting equations 
are, of course, still symmetric with respect to the 
inner product (j^), but split into two sets, com- 
prising the constraint equations and the evolution 
equations, respectively 

2. Among the 10 components of Sg a /3 , 4 correspond 
to infinitesimal coordinate transformations. Hence, 
one has to fix this gauge, or better, to construct 
6 gauge-invariant amplitudes. For a static back- 
ground the components of the linearized extrin- 
sic curvature tensor are gauge-invariant, up to a 
reparametrization of time. If, in addition, space- 
time is spherically symmetric, full gauge-invariant 
quantities can be constructed from the components 
of the extrinsic curvature tensor. 

3. The perturbation equations will be put into the de- 
sired form (0) , where the spatial operator A is first 
required to be elliptic. In order to achieve this, 
ideas introduced by Choquet-Bruhat et. al. M are 
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used. Extending their method, the spatial operator 
is finally made symmetric with respect to the inner 
product (m , where g afl g l3u is replaced by its spatial 
part g lk g^. 



A. The ADM equations 

The Arnowitt-Deser-Misner (ADM) formalism pro- 
vides a 3 + 1 decomposition of (the portion of) space- 
time with topology M = M x S, where S is a three- 
dimensional Riemannian manifold. The manifold M is 
foliated by a one-parameter family of embeddings : 
£ — > M, t 6 M, where the hypersurfaces E t = et(S) are 
assumed to be space-like. With respect to this foliation, 
the timelike vector field d t is decomposed according to 
dt = an + 0, where n is a normal unit vector field or- 
thogonal to S t , and is tangential to Sj. The metric 
assumes the form 

g = -a 2 dt 2 + g l2 \dx l + (3 i dt)(dx j + j dt), (4) 

where the x % are local coordinates on E, and g is the 
induced Riemannian metric on E t . With respect to the 
metric ([|), the Einstein tensor becomes 

G 00 = ^(R-G^ kl K ij K kl ), (5) 
G i0 = G ijkl V j K kl , (6) 

and 

Gu = Gij - 2K\K sj + KKij + ig iS (3K rs K rs - A 2 ) 

- -G. y w V fc Via + -<V'(d t - Lfj)K u , (7) 
a a 

where the index zero refers to the normal vector field 
eo = n = —(dt — 0). The Ay denote the components of 
the extrinsic curvature, 

(8) 

where A' is the trace of Ay, A = g^Kij. All quantities 
with a bar refer to the Riemannian metric g. The tensor 
Qijki j g £jj e p3 e Witt metric on the space of symmetric, 
positive definite matrices, 

Equations (0) and (|J) are a set of evolution equations for 
the 3-dimensional metric g and the extrinsic curvature 
K, whereas eqs. (§) and (§) are constraint equations. 
The freedom to choose the slicing is reflected by the fact 
that there are no evolution equations for the lapse, a, 
and the shift, 0. The Bianchi identities, 



= V^G"" = doG a " + d k G kv + T^ a G a " + T^G^, 

(9) 

guarantee that the Hamiltonian and the momentum con- 
straints, eqs. (||) and (^J) respectively, propagate. 

A convenient choice for the lapse and the shift is given 
by the gauge conditions 

K = const., = 0, 

for which the ADM equations simplify considerably. (For 
asymptotically flat space-times this implies the maximal 
slicing condition A = 0.) In particular, one has for con- 
stant A 

| (fGjj - 3G 00 ) = (A - R - A 2 ) a. (10) 

In vacuum, this yields an elliptic equation for the lapse a, 
which, in a numerical evolution scheme, has to be solved 
after each time step. 

B. Coordinate-invariant quantities 

In a static space-time it is convenient to choose an 
adapted slicing (i.e. orthogonal to dt), such that (3 — 0, 
dta — and dtjjij — 0. As a consequence, the extrinsic 
curvature tensor vanishes, 



The metric then becomes g = —a 2 dt 2 + g, while the 
ADM equations reduce to 

It 1 - - , , 

i?oo = -Aa, Roj = 0, Rij = i?y - -VjV 3 -a. 11 
a a 

Now consider a tensor field Ta V on M, which is static 
on the background, T„ v = 0, T t j = T it = 0. Under an 
infinitesimal coordinate transformation, x^ t— > x 11 + Sx^, 
generated by the vector field A M = 5x^, 6T flv transforms 
according to ST^ ^ ST^ + L X T„ V . Writing (X^) = 
(A 4 , X 1 ) = (f,X% one has 

5T tt ^5T tt + X l T tt . l + 2fT tu 

ST tJ h- 6T tj + X% + fjTtu (12) 
STij STij + X S V S T %] + T is ViX s + T SJ V t X s , 

where here and in the following a dot denotes differ- 
entiation with respect to t, and fj = djf. Using 
Sgu — —2a5a and Sgtj = 60 j, as well as the formulae 
a 2 5T ao = 8T U - 2T tt 5a/a, aST 0j = 5T tj - T ij 80 % , the 
metric and the Einstein tensor transform according to 

6a i— > 8a + aiX 1 + af, 
5f3 3 i-> 60j + Xj - a 2 fj , (13) 
8gij h-> Sgtj + VjAj + VjXi , 
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and 



<5G 0j ^ SG oj + a (Goo/, + G„f ) , (14) 
5G« i * + A -Vj;,., + G ls V 3 X s + G.jViX'. 

The above transformation properties imply that one 
can construct "vector- invariant" quantities, i.e., quan- 
tities which are invariant under the subset of transfor- 
mations generated by vector fields (A M ) — (f,X l ) with 
/ = 0: 



SG, 



SG 



oo — ^Goo — Goo^dp 1 , 

(inv) _ xn 
0j =OLr j, 

sg lj - srvj;,, - G.x.sr - c.x.s.r. 



(inv) 



and 



5a {mv) = Sa- a^Sfij , 



S 9ij 



(15) 



The fact that Sg^J 1 ^ = 2aSKij suggests that - from the 
ADM point of view - the natural gauge-invariant pertur- 
bations on a static background are the extrinsic curva- 
tures, rather than the metric fields. We are, therefore, 
looking for a symmetric wave equation in terms of the 
perturbed extrinsic curvature tensor SKij . 



C. The linearized ADM equation 

The derivation of the linearized held equations is con- 
siderably simplified by choosing an appropriate gauge. 
For a static background the vector field X^ in Eq. (|13| ) 
can be arranged such that Sft = 0. (Note, however, that 
this gauge is not adapted to describe stationary pertur- 
bations.) In this gauge, the perturbed quantities 5a, 5g 
and 5G UU coincide with their vector-invariant counter- 
parts constructed above. The remaining residual gauge 
freedom is generated by the vector field {X^) = (/, X 1 ) 
with Xj = a 2 fj, 



Sa i — ► Sot ~1~ o?oijP + af, 



(16) 



Here and in the following, we use the notation 2oj^ = 

U>ij + U)j i . 

Using the fact that quadratic terms in K give no con- 
tributions if the background is static, the linearized ADM 
constraint equations become 

SG 00 = \SR, 5G l0 ^G ijkl ^5K kl , (17) 

while the evolution equations are 



6Kij = 2a' 9ii ' 



(18) 



SGa = 5G tl - S -G^VfcVza + -G n kl SK kl . (19) 



A wave equation for 5K is obtained as follows: One first 
differentiates Eq. ( [l9|) with respect to t and uses Eq. 
( |l8| ) to eliminate 5g. Applying the general formulae 

25R a0 = - VV M Sg af3 - V a V (g^Sg^) 
+ VV/j 5g au + V^V a Sg Plt 

and 

5G af3 = 5R aP - ^RSg aP - ~g a p (g^SR^ - R^Sg^) 
to the variation of the Einstein tensor G^ in three dimen- 



sions, and also using the background equations (11), then 



yields the following expressions in terms of the vector- 
invariant amplitudes: 



S 



L i:j = aSKij = -Sgij, A = 5a, 



d t S(V=gGi 



(20) 



-kl 



S® G + G t kl UL kl + 2aR%L j)k - 2aR kU] L k 
2V (l aV% )/c - 2X7 k (a (i L j)k ) + 2a k \7 (l L ])k 
2 (V (l V fc a) L j)k - aViVjL - g l:j V k V l (aL k i) 

(21) 



a \a 
where the d'Alembertian is defined according to 

□ = -d t 2 - V r aV r , 
a 

and where 

Sf G = aGij L + gij aG kl L kl - aRLy 
+ GijA. 

In addition, we consider the equation 
-1 



gtjaRL 



-9 



■.dtSlV^gG, 



(HI I 



(22) 



(23) 



-G m A - -G^V^V; ( + G»L 



As expected from the general formula (||) , the spatial part 
of the operators in (Soo, Sy) is symmetric in u = (A, Lij) 
with respect to the inner product 

(«« , u< 2 >) = f + fgfiLtfL®) fj. (24) 
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Furthermore, one has the following constraint equation 
for Li a: 



■ tH 

6G l0 = Gy fei V J ( 

a 



(25) 



induced by the indefinite de Witt metric G % i kl . This pos- 
sibility has been considered in Jll[] . 

In order to get a positive kinematic energy metric, it 
seems unavoidable to choose the maximal (or constant 
mean curvature) slicing condition, 



The system Sy still has the following drawbacks: First, 
the spatial operator is not elliptic: Considering, in partic- 
ular, perturbations of the form L k i = ^Lgij, the second 
derivatives yield — j (VjVj — <7yA)L, which is not gov- 
erned by an elliptic operator. Second, the "kinematic 
energy metric", G' 3 * 1 , is not positive. Indeed, for a sym- 
metric tensor field iy — ty + gPyt, where iy is trace-less, 
one finds Gy- ijy = iy — §<7yi. Finally, there exists no 
evolution equation for the perturbed lapse A. 

In order to overcome these problems, we first apply 
ideas introduced in || , where it was shown that the spa- 
tial operator becomes elliptic after using the momentum 
constraint ( p5| ) and its spatial derivatives. Slightly chang- 
ing the procedure outlined in ||, this can be achieved 
without destroying the symmetry of the final operator. 
It turns out that the combination 



T ■ 

9 ^i. 



0. 



(This gauge may not always exist if the background man- 
ifold £ is closed. In fact, Eq. (|l^) shows that one has 
to solve an elliptic equation for / in order to get L = 0.) 
Adopting the maximal gauge, the trace-less part of Ay- 
yields a manifestly hyperbolic, self-adjoint operator for 
the symmetric, trace-less tensor Ly. However, the per- 
turbations of the lapse are still present in this gauge. 
The amplitude A obeys an elliptic equation which corre- 
sponds to the trace of Ay ■ Hence, the system is of mixed 
hyperbolic-elliptic type. However, as will be shown be- 
low, all terms involving A disappear when the operator 
in Ay is protected onto the constraint manifold defined 
by ©. 



A; 



= d t S {V^gGij) - -V (i (a 3 SG j)0 ) 
-g v 'a 



—g^ k {aHG 



kO I 



(26) 



SP 4 G , where 



almost yields the desired result: Ay = Ay , , , 

A£ = G^ULu + 2aR\ t L ])k - 2aR ka] L kl 

+ 4V (l (a k L j)k ) - 4a {l V k L j)k - 2aV fe (^) L j)k 

- -W i {aa j )L--g ij W k {aa l )L H - 2ag l ^ k {—\ L 
a a \ a / 



1 



G, 



kl 



1 



-.4 



(27) 



(The main difference between Ay' and the f2y introduced 
in H lies in the symmetrizing terms.) 

In order to obtain a symmetric wave equation, one 
makes use of the residual gauge freedom (|l^) and adopts 
one of the following two gauges: The first possibility is 
to choose the harmonic gauge 

A = aL, 

in which case Ay becomes a symmetric operator in Ly . 
Furthermore, this also solves the lapse-problem. Unfor- 
tunately, the De Witt metric still appears in front of the 
d'Alembertian □, and, as a consequence, the spatial op- 
erator is not elliptic. Alternatively, one can multiply Ay 
from the left by the inverse of the De Witt metric, 



r- 1 



1 , 

2 \9ik9jl + : 



'il9jk - 9ij9kl) 



which yields a hyperbolic wave-operator, which is, how- 
ever, only symmetric with respect to the inner product 



D. The linearized Bianchi identities 

In the following the linearized version of the Bianchi 
identities (ji]) will be needed. In particular, we want to 
show explicitely that the constraint equation (p5| ) propa- 
gates. The ADM metric @j with (3 = has the following 
Christoffel symbols: 

r* — — r* — a i r* — — K 

a J a J a 



T k tt =aa k , Y k =aK% T%=t% 



(28) 



Using this in Eq. (j^) yields the 3+1 decomposition of 
the Bianchi identities: 



= V^G 



fj I, 



1. 



= -Goo - aKG oa + -V k (a A G ok ) - aG kl K M , (29) 
a 

= aVG M 

= -G 0j + G ooaj - GojdK + V k {aG kj ). (30) 
Linearizing Eq. ( p9| ) on a static background yields 

SG Q0 - -V k (a 2 5G ok ) = ~aG kl 8K kl ~ aG 00 SK, (31) 
a 

which shows that no new constraint equations for SKij 
are obtained when the Hamilton constraint is differenti- 
ated with respect to t. 

Using Eq. (|3l"|), the linearized version of Eq. ( |30| ) 
becomes 

= OCj - aR jk C k + 2Vj (a k C k ) - 2 aj V k C k (32) 
- aW k (— ) Cj - -W k (a aj )C k - aW k A kj + (BG), 
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where Ay is given in ( f26| ) 
are denned according to 



The constraint variables C, 



SG M = V j 



U 



Cj = a 2 5Goj , 

and (BG) denotes terms which are proportional to com- 
ponents of the Einstein tensor. The fact that Eq. j3^ ) 
also holds for the stress-energy tensor if the matter equa- 
tions are satisfied, implies that these terms need not be 
considered any longer, provided that 8G^ V is replaced by 
<5G>„ — SirGST^ in Eqs. ( |3^ ) and the definition of the 
Cj- 

Hence, if the evolution equations, Ay = 0, and the 
background equations are fulfilled, it follows that the con- 
straints Cj propagate. Note that the operator governing 
the Cj in Eq. (32) is hyperbolic and symmetric. 



The constraint variables, 

C 3 = a 2 (SG oj - 8irGT 0j ) , 

propagate by virtue of Eq. (^) . 

Finally, the differentiated Hamilton constraint yields 



Ann = 



-1 



—Goo A - -V l a 2 V J ) + G»Uj, (35) 
a a \ a 



and the system (Aoo,Ay) is symmetric in (A, Lij) with 
respect to the inner product defined in (|24|). 



E. Summary 

The linearized Einstein equations for a static back- 
ground yield a hyperbolic, formally self-adjoint opera- 
tor for the variation of the extrinsic curvature. This 
quantity is coordinate-invariant with respect to infinites- 
imal diffcomorphisms within the slices orthogonal to dt- 
The wave operator is constructed from the combina- 
tion Aij defined in (26). In the maximal slicing gauge, 
SK = 0, this yields (using also the background equation 
R 0Q = j-Aa) 



klo(BG) 



1V k — 9 &kl 
and 



■2(A-R 



oo 



A V (aa )L k i 

a 



(33) 



A,,- = S. 



(BG) 



i J 

4V 

1 

a 
1 



uLi 



2aR f : L 



(a%)k) 



2aRkUjL 
' a 



4a (l V fc J L, )fe + 2aV fc (^) 



V (l a 2 V 



j) 



;9n 



'-W k {aa l )L kl +AA-RooA) , 



(34) 



where Ay denotes the trace-less part of Ay, and where 
the trace-less symmetric tensor Ly 
defined by 



and the scalar A are 



L, 



t5K 



A = 5c 



These quantities are subject to the residual gauge trans- 
formation (|l|) with V fc (a 2 V fc /) = 0. The system 
(Ay , A k k ) provides a closed set of evolution equations for 
(Lij , A) , which is manifestly hyperbolic and symmetric 
in Lij and elliptic in A. The operator — A + i?oo in equa- 
tion d33| ) is non-negative if the strong energy condition 
and the background equations hold, i.e., if Roo > 0. 



Furthermore, 
equation 



L^ obeys the momentum constraint 



III. THE VACUUM EQUATIONS 

We start by analyzing the wave operator for vacuum 
perturbations of a static vacuum space-time. By virtue 
of the constraint equation 



o = v J I — 

a 



(36) 



the equation Aoo 
rewritten as 



is fulfilled and Eq. (|33|) can be 



AA = 2V k (a l L kl ). 



(37) 



The Laplacian on the left-hand side (LHS) is elliptic and 
symmetric on the dense subspace of C°°-functions in the 
Hilbert space L 2 (T^,fj). As a consequence, its image is 
equal to the orthogonal complement of its kernel. If £ is 
compact, the kernel ker A is the set of constant functions 
on S. If E is not compact, but all perturbations vanish 
sufficiently fast at space-like infinity, ker A is trivial. In 
either case, the RHS of ( |37| ) is orthogonal to ker A and 
therefore lies in the image of A. This shows that equation 
( |37| ) is solvable. A solution, which is unique up to the 
addition of an element in ker A, can be written formally 
as 



A = 2A V k (a l L 



ki ) 



(38) 



where A denotes the restriction of A on the orthogonal 
complement of ker A. 



A. The initial value formulation 



A necessary condition for Einstein's equations to hold 
is that the evolution equations Ay = and the constraint 
equations d36| ) are satisfied. In order to find a solution of 
Einstein's equations, these equations must, however, be 
supplemented with 



G 



= <5G o 



IsR = J-G ijkl S/ ia 2 V-i 



= SGij = 5G l3 



2a 
1 



<5-G,, fei VfcV;a + -^L. 



1 



a 



2 V 



(39) 
(40) 



The first equation is the linearized Hamilton constraint, 
while the second equation is the evolution equation which 
earlier was differentiated with respect to t to construct 
the wave operator. 

The linearized Einstein equations are solved as follows: 

1. Specify any 3-metric Scjij = 5g(t = 0)ij and any 
symmetric, trace-less tensor field Lij on an initial 
time-slice So, such that the Hamilton-constraint 
(p9j) and the momentum constraint (B™ hold. 



2. Choose any convenient function da = 6a(t = 0) on 
So- 

3. Compute A and on Eo from equation (|38|) and 
(ff0|), respectively. By virtue of the background 
equations and the linearization of the Bianchi 
identity (30), Lij/a is trace- and divergence-free. 
Therefore, the momentum constraint propagates as 
a consequence of Eq. ( |32| ) . 

4. Evolve (A, L^) via the equation Ay = and equa- 
tion (|38|). 

5. The linearized metric coefficients are eventually ob- 
tained from 



5a{t) = Sa{0) + J A(r)dT, 
o 

t 

$g(t)ij = S§(0)ij + J 2L(r) 2i dr. 
o 

Note that in contrast to the full non-linear constraint 
equations, the linearized Hamilton and momentum con- 
straint decouple. (This follows from the fact that K van- 
ishes on the background.) The "linearized" conformal 
method to solve the Hamilton constraint is to split 



5ga = 



in its trace and trace-less parts. Inserting this into Eq. 

(H) y ields 



-V fc a 2 V fc | ^ 



2 a \ a 



Again, the RHS of this equation is orthogonal to the 
kernel of the operator on the LHS, which is symmetric 
and elliptic. Therefore, the trace-less part ?fiy can be 
specified freely, and the above equation can be solved for 
</>. How to solve the linearized momentum constraint is 
explained below. 



B. The projection onto the constraint manifold 

Our next aim is to extract the "pure dynamical degrees 
of freedom" by solving the momentum constraint, and to 
project the wave operator onto the space of dynamical 
variables. For vacuum perturbations this can be achieved 
by using the "York decomposition" |[2| . 

First, we replace by 

I 



and the operator A^ by y/ahij . This eliminates the fac- 
tor l/a in front of the second time derivatives in Ay. 
The linearized momentum constraint now reads 



0. 



Let T denote the space of all C°° symmetric, trace-less 
covariant tensor fields on S, and let V denote the space 
of all C°° vector fields on E. We assume that either E 
is compact or all tensor fields vanish sufficiently rapid on 
the "border" of E. Consider the operator 



V, 
2V ( 

la 



and let W € V and Z E T . With respect to the inner 
products on V and T induced by §ij , we have 



(W , L f Z) = (LW, Z), 
■ T is given by 

1 V l W ] +V ] W l -\g lJ S7 k W k 



where L : V 
{LW) l3 -- 



The operator L'L : V — > V is symmetric, positive semi- 
definite and elliptic. Its kernel ker(L^Zz) = ker L consists 
of all conformal Killing vector fields on E. As a conse- 
quence, we can solve the equation 



L^LW = VZ. 



for given Z G T, 



W = (i f i) 1 L t Z- 



Wn 



where (L'L) denotes the restriction of L on the or- 
thogonal complement of its kernel, and where Wq € 
ker L. 

The orthogonal projector on the constraint manifold 



C = (Z i:i ) € T 



= > = ker£ T 



can therefore be represented as 
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P = t-L(L j L) 



Pd?+A)Z 



Any symmetric, trace-less tensor field ZeT can be de- 
composed according to 



7. . - yTT 

4*3 — ^ij 



7TL 



where Z TT /^ 
free, and where 



r,TL 

A i3 



PZ/yfa is trace-less and divergence- 



is the trace- less longitudinal part. This decomposition is 
covariant and orthogonal with respect to the metric on 
T induced by cjij. For a = 1, it corresponds to York's 
decomposition for symmetric tensor fields Kij. York's de- 
composition has further conformal properties which play 
an essential role in solving the full non-linear constraint 
equations [T^ ]. 

In our case, the linearized momentum constraint is 
solved by taking 

Zij P Zij , 

where Zij is any symmetric, trace-less tensor field on 

E. Of course, this requires the computation of (L^L) , 
which is a non-local operator. 



Consider now the wave operator y^aA^ defined in (34). 
It has the form 

yfckij = (d t 2 + A) 

where A is a spatial, elliptic and self-adjoint operator and 
the functions bij (A) are given by 



1 



:Vr,-orV.; 



-9ij 



AA. 



(41) 



By virtue of the linearized Bianchi identity (|32| 
that the spatial operator defined by 

Zij h-> AZij + hj (A), 



it follows 



where A is given by (|3J]), maps the constraint manifold 
C into itself. (In particular, this implies that the mo- 
mentum constraint propagates.) As a consequence, we 
have 

APZ l3 + b %3 (A) = P (APZ tJ + b^ (A)) 



(still provided that equation ( |38[ ) holds). Now it easy to 
see that the bij map into a space which is orthogonal to 
the constraint manifold. Indeed, the spatial operator in 
the system y / a(Aoo, Ay-) can be shown to have the form 

&t 

b A 

As mentioned above, b^P = 0, also implying that Pb — 
0. Thus, the wave equation restricted to the constraint 
manifold takes the form of a symmetric wave equation, 



0. 



where A = PAP is a (formally) self-adjoint spatial op- 
erator. 

In conclusion, we have shown that the wave opera- 
tor, projected onto the (momentum) constraint mani- 
fold, yields a symmetric wave equation for the symmetric, 
trace-less tensor Z+j = i/aSKij. As expected, the varia- 
tion of the lapse does not appear in this equation. The 
following two difficulties remain: First, the quantity Z^ 
is not coordinate-invariant with respect to reparametriza- 
tion of time. Hence, we have not yet isolated all physical 
degrees of freedom. The second difficulty is to find an ex- 
plicit characterization of the constraint manifold in order 
to have an explicit representation of the operator A. We 
have not yet solved these problems for an arbitrary static 
background. However, the additional structure provided 
by a spherically symmetric background, enables one to 
solve both problems. This is discussed in Appendix [A|, 
where a natural derivation of the Regge- Wheeler and Zer- 
illi equations is given. 



C. Locally flat space-times 

To conclude this section, we briefly specialize to the 
case where a = 1. Using the background equations (pi]), 
this implies that (E, g) is locally flat. Equation ( |3^ ) then 
yields A = A(t), and the pulsation equation reduces to 



A ) L 



= Ay = (d t - 

with the constraint 

= V j Lij . 

These equations are well-known from the weak-field limit 
of general relativity. They can be solved by Fourier trans- 
formation. In particular, it follows that all vacuum space- 
times, which can be represented by / x E with I C ]R 
and E locally flat, are linearly stable. 



IV. THE COUPLING TO YANG-MILLS-HIGGS 
FIELDS 

An important new feature of the formulation presented 
in this paper is the fact that it extends to gravitating 
matter fields in a natural way. As mentioned earlier, 
this is due to that - in contrast to the traditional metric 
approach - the wave operator governing the variation of 
the extrinsic curvature appears already off-shell. 

In this section, we show that the equations govern- 
ing fluctuations of a static solution to the Einstein- 
Yang-Mills-Higgs (EYMH) equations assume the form 
of a symmetric wave equation. While this is explicitely 
established for the triplet case, we emphasize that it 
holds true whenever the gauge group is a compact Lie 
group. (Compactness guarantees the existence of an Ad- 
invariant scalar product on the Lie algebra.) 
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A. The ADM equations 

In the ADM formalism, the gauge potential A is 
parametrized in terms of a scalar field $ = — ig n A and a 
one-form A (both Lie algebra valued) , 

A = — $ adt + A, (dx l + /3 l dt) . 

Similarly to the extrinsic curvature, = ^L n gij, in the 
gravitational case, the electric one-form is defined by 

E = -i n F, 

where F = dA + A A A is the YM field strength. We also 
define the momentum belonging to the Higgs field H by 

n = i n DH, 

where D = d + [A\ } denotes the covariant derivative 
with respect to A. 

In terms of these quantities, the ADM decomposi- 
tion of the YMH equations yields the following equations 
(which we list here only for vanishing shift, (3 — 0, since 
they will be needed only for a static background and in 
a gauge with 5(3 — 0; the general case is obtained after 
the substitution dt i— > dt — Lp.): The Gauss constraint is 

-(*D * F - J) = *D*E- [H,U]. 

The evolution equations comprise the definition of E and 

n, 

E = -- (d t A + D(a$)) , 
a 

a 

respectively, and the equations 

- a(*D *F-J) l = -^= d t (VgE 1 ) - *D* (aPf 

- [a<5>,E l ] - [aH,D l H], 

a {*D *DH + 2V'H) = — d t (V|n) - *D* (aDH) 

- [a$,IT] + 2aV'H. 

All quantities with a bar refer to the 3-metric g and to 
the magnetic part of the potential, A. Also, E % = g^Ej 
is raised with the 3-metric. The Higgs potential V is 
assumed to be a function of \H\ 2 = Tr(H 2 ) only, with 
derivative V. 

There is no evolution equation for the electric poten- 
tial <!>, which reflects the invariance of the theory with re- 
spect to gauge transformations of the gauge potential A. 
In fact, $ plays a similar role as the shift f3 for the gravi- 
tational ADM equations. It is easy to see that the Gauss 
constraint propagates as a consequence of the identity 
D 2 * F — [F\ * F] — and the corresponding consistency 
condition D* J = —[DH\ * DH] — [H,D * DH] = for 



the matter current, which holds by virtue of the Higgs 
equations. 

In order to solve the coupled EYMH equations, the 
3+1 decomposition of the stress-energy tensor T^ v — 

T<Z M) +TW is needed: 

+F lk F J k -^g lJ F kl F kl }, 
T ™ } = i Tr U n2 + \(D k H){D k H)} + i-y, 
4 H) = ^Tr{n(Ai?)}, 

T l H) = ^{{DMDjH) ^g~ij(D k H)(D k H) 

Here, Tr denotes an Ad-invariant innerproduct on the 
Lie algebra. The complete set of EYMH equations in the 
ADM formalism is now obtained from the expressions rts), 
(^) and (Q) for the Einstein tensor, and the definition (g) 
of the extrinsic curvature. 



B. Gauge-invariant quantities 

For a static background the slicing may be chosen such 
that g — —a 2 dt 2 + g and Kij — 0. As a consequence, 
= 0, and it is consistent to set E = and n = 
0. Hence, we are considering a static, purely magnetic 
background: 

where A and H do not depend on t. The YMH equations 
reduce to 

*L>* (aF) + [aH, D'if] = 0, 

(aDH) - 2aV'H = 0. (42) 

The perturbed amplitudes, <5<I>, 5E, SH and fll are 
subject to both coordinate and gauge transformations. 
Since E and n vanish on the background, we except them 
to be gauge-invariant to first order: Under infinitesimal 
coordinate transformations generated by a vector field 
X = (f,X l ), one finds, using the same formulae as in 
section [II B| : 

SEi^SEi + aFijf, 
OTh-*. 8IL + a(DjH)f. 







Therefore, 5E and 5H are vector-invariant quantities. 
Under an infinitesimal gauge transformation generated 
by a Lie-algebra valued scalar x, we nn d that 

(5$ M> - — x, 
a 

whereas 5E and 5H remain invariant. Hence, the gauge- 
invariant quantities to evolve are the perturbed electric 
one- form and the perturbed Higgs momentum. 



C. The linearized equations 

In order to simplify the derivation of the wave opera- 
tor, it is convenient to choose the gauge function \ such 
that 5$ = 0. In this gauge, the linearized ADM Gauss 
constraint becomes 

-5(*D * F - J) = *D*5E - [H, 511] , 

while the ADM evolution equations are 

SE = —8t 6 A, 5U = -d t 6H, 

a a 

and 

-aS(*D *F-jy = 5E l - (5*5* (af)' - 5[aH, &H], 

aS (*D *DH + 2V'H) =511- 5*D* (aDH) (43) 
+ 25 (aV'H) , 

These equations have exactly the same structure as the 
gravitational equations (|l7|), ( |l8|) and (|l9|), where 5A it 
5H, as well as 5E i: 5H, 5Kij correspond to each 
other. In order to obtain a wave equation for 5E and 511, 
we therefore differentiate the two equations in (43) with 
respect to t. Again defining vector-invariant quantities 
according to 



£ = a8E = -d t 5 A, # = -a5U = -d t 5H, 
and also using 

d t 5g~ij = 2Lij, d t 5a = A, 
yields the following expressions: 
-ad t 5{*D * F — J)i 



1 



S\ + *D*aD8i - a[F l3 , £ j ] + a [H, [£ h H]] 



aDi [H,V] -2a[DiH,V] 



(44) 



2D j (L kj aF\) + 2aF kj V j L kl - aF i3 V j [L+- 



-ad t 5 (*D *DH + 2V'H) 

- iD*aDt> + iaV"Tr(H^)H + 2aV't> 



[D k (a£ k ),H]-2a[£ k ,D k H] 



(45) 



2D 1 [LijaD^H) + a(DjH) [L+-], 



A 



for the linearized YM equations and the linearized Higgs 
equations, respectively. Here, we have also used the back- 
ground equations (42). Equations (E3) and (EBI) provide 
a system of evolution equations for £ and 'J, which is 
manifestly symmetric with respect to the inner product 



((£«,*«), (£( 2 \*( 2 ))) 

= /Trj^fff+^M 2 ^. 



(46) 



However, we are again faced with the problem that the 
evolution equation governing £ is not hyperbolic. Like 
in the gravitational case, a hyperbolic wave operator can 
be constructed without loosing the symmetry. Again, the 
strategy is to add suitable combinations of the linearized 
Gauss constraint, 



a5(*D * F — J)o = a*D* 



[H,V], (47) 



and its spatial derivatives. Furthermore, one can use the 
linearized momentum constraint, 



d = a 2 6(G M - 8itGT i0 ) 



a 2 G ijkl V j 



(48) 



2GoTr (F k £ k - {D t H)^) 



in order to eliminate the divergence terms of in the 
last lines of Eqs. ( (44| ) and (|4^). The motivation for this 
(besides eliminating first order derivatives) is the follow- 
ing: We will see that the adjoint of the operator governing 
the matter perturbation in the linearized stress-energy 
tensor contains no divergence term. In order to formu- 
late the linearized EYMH equations as a symmetric wave 
equation, no divergences of must therefore appear in 
the YMH equations. The correct combinations turn out 
to be 

A<™> = - a d t 5(*D * F — J)i 

+ -Do?5{*D * F - J) + 2F\C k , 



A (ff) = -adtS {*D *DH + 2V'H) 

- a 2 [5{*D *F- J) ,H] + 2{D k H)C k . 



(49) 



This yields the following hyperbolic, symmetric equa- 
tions: 



A 



(YM) 



and 



□^-a[Fy,^]+a[JJ,[^,if]] 

4GoTr {F l k £i - D k H^) F k - 2[D l {aH), tf] 
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2aP k W k L tJ - -L k3 &{a 2 F kl ) 



+ aF i7 -V J \ L- — \ + 2a k F ki L, 



(50) 



A (ff) = \3V + 4aV"Tr(Hy)H + 2aV'V + a[[H,$>],H] 
- 4GaTr (F l k £ t - D k H^>) D k H - 2[£ k , D k (aH)} 



Lj k D J 



2 D k H) 



A 



a(DjH)V 3 [L + 2a k {D k H)L 



where the d'Alembertians are defined by 



(51) 



1 



1 



-d?£i + D^aDSi + —Da 3 D^ — 



a 



9 t 2 * + D^aDV, 



and where - for Lie-algebra valued p-forms u> on M - 
the operator adjoint of D is D^ui = (—1) p *D*lu. In or- 
der to complete the pulsation equations, we still have to 
compute the linearized stress-energy tensor: Using the 
background equations and 

Sp^ = -DiSj + DjSi, SH = etc., 

one finds that the expression for Ay - defined in (|2 
with G p „ replaced by G M „ — iirGT^ - is given by 

1 



A; 



A 



4GTr < D k aFn £ 



(aF { k £ 0) )+-£ k D { ^F k ) 



1 



£ l a k F kl + -D (l (a 2 D 3) H) * + g ija k (D k H)V 



+ AGaTi ^FIF^Lh 
- 4Ga$L, 



1 



J— 4 
2g ij a^L, 



F u F kl U 



-g %J F kl F kl L 



(52) 



where Ay is defined in Eq. (p7[). The matter perturba- 
tions arising in Ay perfectly fit together with the grav- 
itational perturbations arising in K^ YM ^> and A^ H \ with 
exception of the terms involving (L— — )■ Hence, one 
either has to adopt the harmonic gauge, A = aL, for 
which the operator is symmetric but the kinetic energy 



is not definite (see section II C ) , or the maximal slicing 
condition, L = 0, for which the operator acting on Ly 
is hyperbolic and symmetric, but the perturbed lapse is 
still present. As in the gravitational case, we adopt the 
maximal slicing condition. 



D. Summary 

The linearized EYMH equations (with arbitrary com- 
pact gauge group) for a static, purely magnetic back- 
ground, yield a hyperbolic, formally self-adjoint operator 



for the variation of the extrinsic curvature, the electric 
one-form and the Higgs momentum. All these quanti- 
ties are invariant with respect to both infinitesimal gauge 
transformations of the gauge fields and infinitesimal co- 
ordinate transformations within the slices Ef. The wave 
operator is constructed from appropriate combinations 
of the linearized ADM equations, as defined in Eqs. fl26| ) 
and (^]). In the maximal slicing gauge, the evolution 
equations are described by the trace-less part of the ten- 
sor Ay, given in Eq. (152), as well as the expressions 



for A< yM ) and A^ defined in Eqs. © and ©, re- 
spectively, where one sets L = 0. The constraint equa- 
tions are the momentum constraint, (^8|) and the Gauss 
constraint, (^7|). (Additional constraints involving also 
perturbations of the metric, the gauge potential and the 
Higgs field themselves are the Hamilton constraint and 
all evolution equations, which were differentiated with 
respect to time in order to construct the wave operator.) 
Furthermore, the trace of the tensor Ay yields the fol- 
lowing elliptic equation for A: 

(A - Rao) A = 2V k {a l L kl + GTv{aF lk £ 1 )} 

+2GaTr { [H, D k H]£ k + 2V ' - F km F l m L kl ) , 

where the momentum constraint and the background 
equations have been used as well. At least in the pure 
EYM case, this equation is always solvable, since either 
2i?oo = GTv(F k iF kl ) > and the operator on the left is 
invertible or the magnetic field vanishes and the equation 
reduces to its vacuum counterpart (|37|). In the presence 
of Higgs fields further investigations are needed, since 



(H) 
00 



-2GV is negative in that case (i.e. the strong 



energy condition does not hold for the Higgs part of the 
stress-energy tensor). In some relevant situations, how- 
ever, the trace equation causes no problems: In particu- 
lar, this equation is void for odd-parity perturbations of 
a spherically symmetric background. 

Finally, it is also instructive to compute the differ- 
entiated Hamilton constraint, defined as in (|35|)- with 
G MI , replaced by G^ v — 8ttGT^. This yields aAoo = 
— V J Cj, where Cj are the constraint variables defined 
in Eq. (B3). The system of equations defined by 
(A o, Ay, Af M \ A( H ^) is symmetric in the amplitudes 
u = (A, Lij,£i, \&) with respect to the inner product 



g lk g 



L i] L kl ' 



+ 2GTr \g^£\ l) £f ) + iJrC 1 )^ 2 )}] fj. (53) 

As Aoo = holds by virtue of the momentum constraint, 
this again implies that all terms involving the ampli- 
tude A in the remaining equations Ay, a[ YM ^ and A^ H ^ 
are orthogonal to the (momentum-) constraint manifold. 
Therefore, the variation of the lapse does not appear if 
the evolution equations are projected on the constraint 
manifold. The initial value p roblem is solved similarly to 
the vacuum case, see section III A. 
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V. THE COUPLING TO PERFECT FLUIDS 



B. Coordinate-invariant quantities 



At a first glance, it may be surprising that our formu- 
lation also applies to perfect fluids. It is, however, known 
pl| that Lagrangian formulations exist in this case. As 
we show in Appendix the linear fluctuations of a self- 
gravitating perfect fluid are described by the linearized 
metric and the Lagrangian displacement vector. If the 
background is static, it will be shown that the wave op- 
erator acts on the components of the extrinsic curvature 
and the time-derivative of the displacement vector. 



For a static background one has G^o = 0, and it is 
consistent to set ui — (i.e. the flow is static). The 
ADM equations then reduce to 



= -Aa - 4 7 rG( / 9 + 3P), 
a 

= Rij - iViVjQ - 4nGgij(p - P), 



o 



= VjP +(p + P) 



a j 



(58) 
(59) 
(60) 



A. The ADM equations 



Under an infinitesimal coordinate transformation gen- 
erated by (X^) = (/, X 1 ) one has the following behavior: 



The stress-energy tensor for a perfect fluid is given by 



= (p + P^^Uu + Pg 



where u is a time-like vector field, describing the fluid's 
motion, normalized such that g(u,u) = — 1. The fields p 
and P denote the energy density and the pressure of the 
fluid, respectively. The ADM split of u is given by 



7(e + v*ei), 



(54) 



where 7~ 2 = 1 — g(v, v) and eo is the previously intro- 
duced normal unit vector field orthogonal to £(. We 
choose the sign of 7 such that the fluid flows in the same 
direction as eo- 

The relativistic Euler equations are obtained from 



jLV 







0: 



IV 



'V7 U 



= q a u V u T av 



(P- 

(P- 



p)e + 



f <z q m v q p, 



(55) 
(56) 



where 



is the expansion, 



(V„u) A 



«"V l; « (1 is the acceleration vector, and where q a p = 
S a p + u a U/3 projects onto the spaces orthogonal to u. In- 
serting the expression ([54J) above yields 

9 = -7 + jK + - V 3 (au, ) , 
a a 

a — (Ui + 70j) + u J VjUi. 

Here, we have assumed that the shift vector (3 vanishes, 
such that Ui = 71;*. For linear perturbations of a static 
background there is no loss of generality in doing so since 
we will choose a gauge with 8(3 = 0. The stress-energy 
tensor gives 



r 00 = 7 2 G°+n 

Tio = -j(p + P)ui, 

= (p + P)uiUj + Pgij. 



(57) 



In order to get a closed system of equations, an equation 
of state is needed as well. In the following we will only 
assume that this has the form P = P(p). 



Sui i * Sui - afi, 
8p^ 8p + X s V s p, 
SP^ SP + X S V S P. 

Therefore, Sui is a vector-invariant quantity, and using 
(|l3|), we may also construct the vector-invariant ampli- 
tudes 

S p (inv) =s p - S(3 s V sP , 

S p(inv) EE $p_ SP S V S P. 

Using the transformation properties of Sui it is tempting 
to construct full coordinate- invariant amplitudes. These 
correspond to Lagrangian deformations of g, p and P 
(see Appe ndix |b[). Moreover, it turns out that using the 
relations (B5), one can find a wave equation for Ag^- 
alone. Unfortunately, the resulting pulsation equation is 
not symmetric. 

The correct quantities to evolve - from the philosophy 
adopted in this paper - are the vector-invariant quantities 
SKij and Su^. This is motivated by the fa ct t hat for 
a static background, the second relation in (B5) yields 
aSui = 



C. The linearized equations 

We start with the linearization of the relativistic Euler 
equations ( ^5|) and (^). As before, we make use of the 
gauge freedom in order to set 8(3 — 0. Noting that 7 is 
quadratic in m and that Ui vanishes on the background, 
we obtain 

Se = SK+-Vi{a8u ), 
a 



8a,i = —Siii 
a 



Also using the background equation (^), the lineariza- 
tion of equation (Ka) gives 



Sp= [a z (p + P)SK 

a 1 



■ V 



k [a 2 (p + P)6u k ]). (61) 
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Defining the speed of sound c s by c 2 — dP/dp, we also 
have 

2 

SP = ~ (« 2 (P + P)SK + V fc [a 2 (p + P)5u k ] J . (62) 



[The above equations also follow from the relations ( JB5| ) 
derived in Appendix Finally, using the fact that P 
is a function of p only, the linearization of equation ( |5C " 
yields 



-Jut 



<5a 



SP 
P + P 



(63) 



An evolution equation for Sui is obtained after differ- 
entiating this equation with respect to t and eliminat- 
ing SP using equation (|62]). Defining £> = a 3 (p + P), 

Wj = a 2 {p + P)Suj (64) 
and recalling that L — aSK , ^4 = 5d, we find 

c^VT^) = ^d t 2 Wi - V, (GV fc W^ fe ) 

+ %( ^ - c*lJ = 0. (65) 

Furthermore, equation (^33|) yields the following circular- 
ity condition: 




- V, 



B 



= 0. 



(66) 



This condition would also allow to introduce a poten- 
tial $, defined by Wi — BVj$. However, in order to 
obtain symmetric equations, it turns out to be neces- 
sary to formulate the pulsatio n eq uations in terms of Sui 
rather than $. The condition (66) is, therefore, kept as a 
constraint equation for the system. That this constraint 
propagates is obvious from Eq. (|6^). The circularity 
condition is used to cast the evolution equation into a 
hyperbolic and symmetric form. Indeed, taking a time 
derivative and using Eq. (p6|), one can rewrite Eq. j65| ) 
in a manifestly hyperbolic form without destroying the 
symmetry: 



= \;d 2 W, - V k CV k W t + F k V l W k - V fc 



B 
1 



FiW k 



B B z 



CR lk W k 




(67) 



where F k = V k (BC)/B. 

In order to complete the pulsation equations, we have 
to compute the linearized stress-energy tensor. Using 
@ and © in the expression for Ay - defined in (Eq) 
with G p „ replaced by — 8wGTfj, v - we find 



A; 



A) 



8ttG 



(a(p-P)L t 



+ -G l]kl V k (aW l ) + (1 - c 2 s )g l3 V k W k 



agij [c 2 s (p + P)-(p-P)]L 



(68) 



where A^ is defined in (|27|). In order to obtain symmetric 
equations, we use the linearized momentum constraint, 



r k 

Ci = a z S(G l0 - 8ttGT i0 ) = afGyuV 3 ( — 



SirGWi. 
(69) 



Defining 

A( F ) = d t S{q^T av ) 
we finally obtain 



-Ci 



A (F) = ~0 t a Wi - V, (CV k W k ) 



16ttG, 



-Wi 



a 



a 



2aGy fe ,V J ( — ) + V, ( (1 - ct)L - L 



This shows that, with the exception of the last term in- 
volving (L — — ), the matter perturbations arising in Ay 
perfectly fit together with the gravitational perturbations 
arising m 

Like in the vacuum case, we now adopt the maximal 
slicing gauge L = 0. The evolution equations for Ly 
and Wi are given by the trace-less part of Ay and A^ F ) , 
where one sets L = 0. The constraint equations com- 
prise the linearized momentum constraint, (|69|), and the 
circularity condition, (|6^). The trace part of Ay yields 
the following equation for A: 

(A - R Q0 ) A = 2V k {a l L kl ) - 4ttG(1 + 3^)V fe W fe , 

with Rqq = AttG(p + 3P). Here, the momentum con- 
straint has been used in order to simplify the equation. 
As long as p + 3P is positive, the operator on the LHS 
is invertible and thus the equation for A is solvable. The 
differentiated Hamilton constraint - defined in Eq. ( |35| ) 
with G M „ replaced by G M „ — iirGT^ - yields 



An 



1 



a 2 V^ [ 

a 



i-n-GWi 



The system of equations defined by (Ann, Ay , A^ F )) is 
symmetric in the amplitudes u = [A, Zy,Wi) with re- 
spect to the inner product 



AWaW + <f k g^L^L^ + SttG^W^W^ 



As before, Aoo = is satisfied by virtue of the momentum 
constraint and, as a consequence, the variation of the 
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lapse does not appear if the wave operator is projected 
onto the constraint manifold. 

Finally, the evolution equations Ay = and 
can be written in a manifestly hyperbolic form by a 
further differentiation with respect to t, as we have shown 
above. The initial value problem is solved similarly to the 
vacuum case. 



VI. CONCLUSION 

We have shown that the perturbation equations gov- 
erning linear fluctuations on a static background can be 
cast into the form of a constraint symmetric wave equa- 
tion for gauge-invariant quantities. In particular, we have 
discussed the initial value formulation for vacuum fluctu- 
ations of vacuum space-times, where the adapted gauge- 
invariant quantities are the components of the linearized 
extrinsic curvature tensor. For a spherically symmetric 
background, the constraint equations have been elimi- 
nated in a natural way, and the equations of Regge- 
Wheeler and Zerilli were rederived. An important new 
feature of the curvature-based perturbation formalism 
presented in this paper is that it admits a natural gener- 
alization to gravitating matter fields. 

As a first example, the pulsation operator governing 
fluctuations on a static, purely magnetic EYMH con- 
figuration were derived for an arbitrary compact gauge 
group. These pulsation equations are expected to be 
valuable in order to discuss the stability of spherically 
symmetric solutions to the EYMH equations with re- 
spect to non-spherical symmetric perturbations. In a 
forthcoming article [|13| , we will show that the Bartnik- 
McKinnon solitons |14f and the corresponding black holes 



with hair 
sector wit 



If 



admit no unstable modes in the odd-parity 
r total angular momentum £ > 1. It should 
also be interesting to generalize the investigations to so- 
lutions with Higgs fields or with a negative cosmological 
constant Jig] , since some of these solutions are known 
to be linearly stable with respect to radial perturbations. 
Also, the stability of new static, axially symmetric config- 
urations |l^] can be discussed within the new framework. 

In a second example, we have shown that a symmetric 
formulation of the equations governing linear fluctuation 
of a self-gravitating perfect fluid exists. Provided that 
the constraint and the dynamical variables can be decou- 
pled in a "symmetric way", the resulting equations are 
expected to be useful in order to study the stability of 
static, relativistic stars, or the emission of gravitational 
waves from pulsating neutron stars |jl8), |Q. In order 
to discuss the fluctuations of rapidly rotating configura- 
tions, it is our aim to generalize the curvature-based ap- 
proach to general stationary, i.e. non-static, background 
configurations. However, this requires further investiga- 
tions, since the derivation of the wave operator must be 
generalized to cases where the extrinsic curvature of the 
background no longer vanishes. 



Finally, we mention that the symmetric formulation 
could also be useful for second (or higher) order perturba- 
tion theory, since the n'th order perturbation is governed 
by the same wave operator as the first oder perturbation, 
but with a source term depending on perturbations up 
to order n — 1. 
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APPENDIX A: DERIVATION OF THE 
REGGE- WHEELER AND ZERILLI EQUATIONS 

In this Appendix, we specialize the general formalism 
to a static spherically symmetric background and show 
how to project the wave operator onto the momentum 
constraint manifold. This will result in a new - and 
transparent - derivation of the Regge- Wheeler and the 
Zerilli equations. 

A convenient parametrization of the spherically sym- 
metric 3-metric is 



dx 2 + r 2 g, 



(Al) 



where a; is a radial coordinate and where g — dfl 2 is the 
standard volume element on S 2 . The field r and the lapse 
a are functions of x only. The Christoffel symbols with 
respect to g become 



0, T X AB = -rr'g A B, 



r xx — o, r lB - —S 



Bi AB 



r° — r 1 

1 A R — 1 



AB ' 



(A2) 



where capital indices refer to coordinates on the 2-sphere, 
and a prime denotes differentiation with respect to x. 
The components of the Riemann tensor are given by 

R X AxB = - rr "9AB, R X CAB = °> 

R D CAB = 2(1 - r' 2 )5 D [A g B]c . 
Einstein's background equations ( |Tl|) become 
(r 2 a')' 



— R 00 

i 

= Rxx = — 2- 



= R A b = | 1 - r' 2 - rr" - rr' — J w . 



(A3) 



The solutions to these equations are, of course, the 
Schwarzschild solutions 



of = 1 - 



2Gr 



dr 2 



r 2 dn 2 , 



where m is constant. Our aim is to write t he w ave opera- 
tor Aij in terms of the background metric (Al). Since the 
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background is spherically symmetric and invariant under 
parity reflection, x i— > — x, we can perform a multipole 
decomposition in the odd- and the even-parity sectors 
separately. As this decomposition can be done in an or- 
thonormal manner, the resulting wave equations will be 
symmetric again. 



1. The odd-parity sector 

The symmetric tensor field can be expanded in 
terms of spherical tensor harmonics. In Appendix D of 
Ref. H we have shown how these tensor harmonics can 
be obtained from the standard scalar spherical harmon- 
ics Y im . Since the background is spherically symmetric, 
perturbations belonging to different I and m decouple; in 
the following we shall therefore suppress these indices. 

In the odd-parity sector, the expansion is given by 



L X x — 0, 

LxB = ni m Sb, 

Lab = n 2 ru 2 2V( A S B ) 



(A4) 



where Sb = {*dY)s denote the transverse spherical vec- 
tor harmonics. The functions u\ and u 2 depend on x 
only, and the normalization constants n\ and n 2 are cho- 
sen such that 



(Lij , Lij) = j (u\ + u\)dx, 


where (. , .) denotes the inner product (p4|). One finds 

m = [2l{i + l)f 1/2 , n 2 = [2i{l + l)Af 1/2 , 

where we have defined A = (£ — l)(t + 2). Since the 
perturbation of all scalar quantities vanishes in the odd- 
parity case, u\ and u 2 are coordinate-invariant and the 
perturbation of the lapse vanishes. As a consequence, we 
obtain a symmetric wave equation for the amplitudes ui 
and u 2 . Introducing ( A4) into the wave equation Ay- = 
defined as in (|34|) gives 



to 



? - 



S + 2V BG ) u = 0, 



(A5) 



where 



L (1) 

7 \ r ) pp 

2%/A 7p 



A 7 2 2VX 



7p 
A 7 2 



Here, we have defined 7 = a/r and the new radial 
coordinate p according to dx — adp. Furthermore, 
Vbg = 7 2 (l — r ' 2 ~ rr " ~ rr'a' la) vanishes by virtue 
of the background equations ( |A3| ) . The momentum con- 
straint equation yields = Cb = tiiu c Sb, where 



7 r 



v/A 



7 "2 



(A6) 



parametrizes the constraint variable. 

Our aim is to find a new variable u p , representing the 
dynamical variable, such that the wave equation assumes 
the form 



d: 



v c 
Vpc v p 



0. 



(A7) 



On the constraint manifold, u c — 0, one then has 

[d t 2 - d 2 p + v p ] Up = 0, 

which is a symmetric wave equation for the dynamical 
variable u p 
the ansatz 



variable u p . It turns out that this can be achieved with 



B 



Ui 
»2 



B = d a 



A 



where the first row of the matrix has been chosen such 



that Eq. (A6) holds, while the second row has been cho- 
sen such that the matrix is symmetric. This guarantees 
that no first order derivatives appear in B. The key 
observation is that the function A can be chosen such 
that 

-d 2 + S = B^B. 
This is indeed the case if A = — — . As a consequence, 



the desired wave equation ( A7) is obtained aft er a pplying 
B to the left of the original wave equation (|A5|). This 
yields 



d 2 



BB^ 



0. 



where 



BB' 




+ 7 2 A 



(?) +7 2 A 

V T I p p 1 



Thus, the constraint variables decouple from the dynam- 
ical variables, and the wave equation governing the dy- 
namical degree of freedom is exactly the Regge- Wheeler 
equation |8|] . It is also worthwhile noting that the spatial 
part of the wave operator can be factorized as B^ B, with 
B regular. This implies that the Schwarzschild solution 
is linearly stable in the odd-parity sector. It turns out 
that this factorization can be generalized in the presence 
of a SU(2) Yang-Mills field. This then implies the ab- 
sence of non-spherically symmetric unstable odd-parity 
modes for the Bartnik-McKinnon solitons |14| and the 
corresponding black holes with hair |15[ ]13[. 



2. The even-parity sector 

In the even-parity sector is expanded according to 
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L xx =ni-Y, 
r 

L x b = n 2 qV B Y, 

T t (3) i (4) 

Lab = -n 3 h e\ B + n 4 5 e\' B , 



(A8) 



where 

= rg AB Y, e% = r ( V , V ;; V + ^(/ 4 1 )/, U ;>" 



The basis is chosen such that Lij is trace-less and or- 
thonormal with respect to the inner product 



(Lij , Lij) 



(h 2 



g 2 )dx. 



One finds 



3, - 2 = 



n 3 = -^, 



2 



where we have defined fj 2 = ^(£+1) and A = (l-l)(£+2). 

In contrast to the odd-parity case, the amplitudes h, q 
and g are subject to residual coordinate transformations 
of the form (|l6|). Expanding / — £(x)Y, we find 



n 3 h^n 3 h~ a r £ + — £ , 

2 r 

n 2 q h-> n 2 (? + ar f— j 



(A9) 



r 

where / obeys the equation V fc (a 2 Vfc/) = 0, i.e. 
1 ■ - 2 



0. 



As the dynamical variables must be coordinate-invariant, 
these residual transformations will help us to construct 
them as approp riat e combinations of the amplitudes 
(h,g,q). Using (A4) in Ay = yields the wave equa- 
tion 



(d 2 -8 2 + S)u + b{a) = 0, 
where the symmetric matrix S is given by 



S = 



with 




12m7p 
12 M7p 7 2 S 22 2 N /A 7p 
2y/X lp j 2 S : 



(A10) 



(All) 



3.3 



S11 


= M 2 


+ 6r' 2 


c " 

— orr 


— 9rr' 


0' 










a 


S22 


2 

= A 4 


+ 4r' 2 


— 4rr" 


— 8rr' 


a' 










a 


S33 


= M 2 


-2r' 2 


- rr" - 


, a 
- rr — 





CV" 



-{aa')' 



The inhomogeneous term b(a), where a = a(x) is the 
scalar amplitude parametrizing A, A = a aY, is found to 
be 



n 3 61 = 2a r'a' — \i — a 
r 

n% b 2 = -2ar (—a) , 

n± 03 = a. 

r 



(A12) 



Next, the Adequation yields 



(r 2 (cra)') — a\j?a 



4^3 
a 



aa' (act)' 

r 



0. 



This equation has already been used in our derivation 
in order to eliminate the second derivatives of a in the 
expression for b\ . Finally, the linearized momentum con- 
straint yields the following two equations: 



= C x = 2n 3 
= C B = n 2 



air 



a j r q 



2 



-7 9 



Y 

r ' 



7 h — VA7 .1 



(A13) 



VbY. 



Unfortunately, we did not succeed in finding a trans- 
formation from the original variables u to the dynamical 
variables u p in a similarly elegant way like in the odd- 
parity sector, i.e. by means of a supersymmetry trans- 
formation. However, the dynamical degrees of freedom 
are found systematically as follows: First, one has three 
original variables u — (h, q, g) and two constraint vari- 
ables, which are given by Eqs. (A13). We are, therefore, 
looking for a single dynamical scalar amplitude. This 
amplitude must be inva rian t with respect to the residual 
gauge transformation (A9). The simplest way to con- 
struct a completely coordinate-invariant amplitude is to 
take a linear combination of the amplitudes h, q and g. 
Requiring invariance with respect to the residual trans- 
formation ( |A9| ), the only possibility (up to rescaling) is 

Z = n 3 h + n 2 r'q - n 4 fg, 

where the function / = f(x) is given by 

f(x) = ^ 2 +rr'--r' 2 . 
2 a 

Taking the corr esponding combination of the evolution 
equations (A.10) yields 



(d t 2 - d 2 ) Z + j 2 [i 2 Z - 6aa'—Z ~ 6n 2 aa'—q = 0, 



where the background equations (A3) have been used 
in order to simplify the expression. The coordinate- 
invariant amplitude q is given by 



1G 



n± ( g 



(-) 



n 2 \a/ P 

Note that the lapse amplitude a does not appear in the 
above equation. This becomes clear when one compares 
( |A12| ) with (|A9|): Any gauge- invariant combination of 
the evolution equations annihilates the terms involving 
a. The constraint equations, (A13) yield 



n 2 fq = 0. 



The amplitude q can therefore be eliminated. After using 
the background expressions 

2 ,o AI „ 2Gm ,r' GmN 
a = r = N = \ 



act — 

r 



we obtain the scalar equation 



A + 



6GV 



12Grn 



N 



(9 t 2 - a p 2 ) Z - —^NZ P + f\-Z = 



Defining Z = (A + 6Gm/r)f , this is the Zerilli equation 

i-. 

Finally, we mention that the separation between the 
constraint and dynamical variables can also be achieved 
by a first order transformation of the form 'J = Ah + 
Bq + Cg + g p with suitable functions A, B and C . This 
ansatz results in the Regge- Wheeler equation pq| . A 
similar result has been obtained in ||, however, in con- 
trast to their derivation, the procedure adopted here is 
more natural, since the requirement of gauge invariance 
with respect to reparametrization of the time coordinate 
implies that only one function must be matched. 



APPENDIX B: PERFECT FLUIDS AS FIELD 
THEORY 

In this appendix, we briefly review a field theoretical 
formulation of perfect fluids models |2l| and discuss some 
applications to linear perturbation theory. 

The fluid's motion is described by a map 

F:M^T,x^(F 1 (x),F 2 (x),F 3 (x)), 

where (M, g) denotes space-time and where T is a three- 
dimensional manifold. The fields Fx, F 2 and F3 are scalar 
fields on M and describe a material coordinate system for 
the fluid. We assume that the linear map F* x : T X M — > 
Tfi x )T has rank three and that the kernel of F* x is time- 
like for all points x <E M. 

Let f2 = n(F)dFi A dF 2 A dF 3 be a volume form on 
r, where n{F) describes the material density. Then, we 
define a 4-current J by 



J 



where F* is the pull-back and * the Hodge-dual on 
(M, g). In local coordinates, we have 



J„ = n( J F) W7 (V Q J F 1 1 )(V 3 J F 1 2 )(V^F 3 ). 

[As an example, consider F a (t, x a ) — x a —tv a , v a — const, 
which describes a fluid in constant motion in flat space. 
Then we have J* = h, J % = nv 1 .] The definition of J 
implies that 

*d * J = 0, 

thus the particle number is conserved. Note also that 
J spans the kernel of F* . The particle density and the 
4-velocity, as measured in the coordinate system x, are 
defined by 

1 



n 

respectively. 

The dynamics of the fluid is described by a Lagrangian 
of the form 



C(F a ,VF a ) = -nu(v). 



(Bl) 



Here, u is a function of v = 1/n and describes the inner 
energy per particle. Some calculations show that the 

> <5£ 



stress-energy tensor, T^ v — + g^C, is given by 



(p + P)u»u v + g^P, 



where p = nu is the energy density and 

du 
dv 



P 



is the pressure. Note that this agrees with the first law 
of thermodynamics, 5u = —PSv, when the entropy per 



particle is constant. Using T^ u 
one also finds 



V„T^ = (d p F a ) 



dC 

dF a 



dC 



d{d v F a ) 



Hence, the relativistic Euler equations, V„T^ = 0, are 
equ ival ent to the Euler-Lagrange equations with respect 
to (Bl), since F* was required to have full rank. 

Consider linear perturbations of a self-gravitating 
fluid, described by the (Eulerian) perturbations 5g pu and 
5R a . One can also introduce the Lagrangian perturba- 
tion operator A = S + L^, where the Lagrangian dis- 
placement field £ is a first order quantity defined such 
that the material coordinates do not change under the 
perturbations: 

AF a = 0, a =1,2, 3. (B2) 
Since F a are scalar fields, this means that 



5F a = -(d^F*)? 



(B3) 



Therefore, £ M is uniquely defined up to a vector propor- 



tional to 



Furthermore, under an infinitesimal coordi- 



nate transformation generated by a vector field A M , we 
must have 
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(B4) 



(modulo a vec tor proportional to u^) in order to mai ntai n 
equation (B3). As an important consequence of (B4), 
Lagrangian perturbations of a tensor quantity are auto- 
matically coordinate invariant to linear order. 

Since AR a = and [d, A] = 0, A J M and therefore also 
An, Ay, 1 *, Ap and AP can be expressed in terms of Ag 
only. Using 



1 



AJ M = --( 5 qP A 5q/3 )J m 



J a Ag a , 



the following relations can be derived in quite an efficient 
way 



An 
Au^ 

Ap 
AP 



1 



nq af3 Ag a 



/3' 



-(u a u Ag a p)u^, 



(p + P)q af3 Ag a 



/3' 



(B5) 



- 7 Pq^Ag a 



Here, we have also defined the adiabatic index 
d log P n 3 d 2 u 



7 



<91ogn P dv 2 



and q a/3 = g a P + u a u^ . These relations, which are well- 
known (see, e.g., p^]), imply that the Lagrangian per- 
turbations of the particle conservation equation and the 
continuity equation are fulfilled. More precisely, we have 



A , , / An 1 

V u n + n6 = Vj + -q ap Ag a0 

n V n 2 



A (V u p+(p + P)0) = V./ A/ ' 
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